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Due to the particular geometry of the kagome lattice, it is shown that antisymmetric 
Dzyaloshinsky-Moriya interactions are allowed and induce magnetic ordering. The symmetry of 
the obtained low temperature magnetic phases are studied through mean field approximation and 
classical Monte Carlo simulations. A phase diagram relating the geometry of the interaction and 
the ordering temperature has been derived. The order of magnitude of the anisotropies due to 
Dzyaloshinsky-Moriya interactions are more important than in non-frustrated magnets, which en- 
hances its appearance in real systems. Application to the jarosites compounds is proposed. In 
particular, the low temperature behaviors of the Fe and Cr-based jarosites are correctly described 
by this model. 



I. INTRODUCTION 

In the last few years, geometrically frustrated antifer- 
romagnets have been the subject of much experimental 
and theoretical works Q|. Up to now, the most exten- 
sively studied are the kagome (D = 2) and the pyrochlore 
(D = 3) antiferromagnets. Both are expected to have 
disordered classical and quantum ground states and be- 
have as spin liquids ||, ||. In such cases, it is expected 
that any small perturbation may have a strong effect on 
the ground sate manifold. These perturbations can arise 
from thermal or quantum fluctuations, anisotropy, longer 
range interactions, etc. This paper focus on the kagome 
lattice where it is shown that Dzyaloshinsky-Moriya in- 
teractions (DMI) may be present. As a consequence, 
magnetic ordering can occur at low temperature. In par- 
ticular the magnetic structures of the Fe and Cr jarosites 
H O, 12, |l3| are explored and it is proposed that DMI 
can explain the low temperature behaviors of these com- 
pounds. 

DMI may be present in magnetic systems if there is no 
inversion center between two magnetic sites. Such inter- 
actions between two sites i and j are defined by a vector 
Djj-: Hij — Dj_y.(Sj x Sj). It was Moriya Q that clarified 
the conditions for the existence of these interactions and 
he gave some rules determining the possible directions 
of Djj . He proposed a microscopic derivation of these 
interactions based on Anderson's formalism of superex- 
change, including spin-orbit coupling. Other mecanisms 
were proposed for metallic systems with RKKY exchange 
interactions || ^| . 

In fact the most general spin hamiltonian for two 
neighboring spin-1/2 magnetic ions is given by: 



Hi. 



JijSi.Sj 



(1) 



where the sec ond term is the antisymmetric DMI and 
the last term Ay- is an anisotropic symmetric exchange 
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interaction. Only antiferromagnetic isotropic exchange 
(Jy>0) will be considered in the following. 

In antiferromagnetic oxydes such as a-Fe203 [Q, 
La 2 Cu0 4 @, DMI is responsible for weak ferromag- 
netism. In some Fe and Cr jarosites, a small ferromag- 
netic component was observed even when the main ex- 
change interaction is antiferromagnetic (Ocw ~-700 K 
for Fe jarosites). The low temperature magnetic struc- 
ture is a long range ordered state where all the spins have 
the same component in the direction perpendicular to the 
kagome plane, giving rise to weak ferromagnetism. The 
in-plane components of the spins form a q=0 structure, 
the three spins of the triangular magnetic unit cell being 
at 120° one from another. Only one of the two possible 
chiralities is observed. Depending on the jarosite, the 
out-of-plane component may vanish || . 

Each of the magnetic (Fe and Cr) atoms which form the 
kagome lattice is surrounded by an octahedron of oxygen 
atoms, and two neighbouring octahedra share one oxy- 
gen atom which mediates the superexchange interaction 
between the magnetic sites. These octahedra are respon- 
sible for the crystalline electrical field on the magnetic 
atoms. They are slightly distorted and their local axial 
axes are tilted with respect to the normal of the kagome 
plane. 

In section |l[ some aspects of DMI specific to the 
kagome lattice are discussed and a microscopic deriva- 
tion of the DMI is made. Sections III and |y| deal with 
the magnetic properties due to the two types of possible 
DMI. The results are compared to the magnetic struc- 
tures of the jarosites in section |v|. 



II. DM INTERACTIONS IN THE KAGOME 
LATTICE 



As a consequence of the hamiltonian being invariant 
under the symmetry operations which leave the lattice in- 
variant, the direction of is geometrically constrained 
and follows rules explicited by Moriya 0. In the next 
section, Moriya's rules arc applied to th e kag ome lattice 
and to the related jarosites. In section II B, an estima- 



tion of the order of magnitude of the different terms in 



2 



([j]) is made. Section II C gives a microscopic derivation 
of DMI following Moriya's formalism and taking into ac- 
count the peculiar structure and the environment of mag- 
netic atoms in the jarosites. 



A. Application of symmetry rules 



Two of Moriya's rules give useful informations about 
the Dy in the kagome lattice. First, the middle point be- 
tween two sites is not a center of inversion for the kagome 
lattice, so DMI are not forbidden by the symmetry of the 
lattice. Furthermore, in a perfect kagome lattice, the Dy 
can only be perpendicular to the kagome plane since this 
plane is a mirror plane. These symmetry considerations 
determine the axis of all the Dy vectors (if DMI exist), 
but not their directions nor their values, which will de- 
pend on microscopic details. 

In the jarosites, the symmetry is lowered because the 
octahedra of oxygen atoms which surround the magnetic 
sites are tilted 12, The local axial axes of these 

tetrahedra are not exactly perpendicular to the kagome 
plane, and the kagome plane is then no longer a mir- 
ror plane for the lattice when the non-magnetic atoms 
are considered. These non-magnetic oxygen atoms that 
make up the coordination octahedra must be taken into 
account because they are responsible for the crystalline 
electric field on the magnetic atoms, and are involved in 
the superexchange mechanism between these magnetic 
sites. Applying Moriya's rules to the jarosites crystal 
constrains the Dy vectors to be in the plane perpendic- 
ular to the bond (ij) since this is a mirror plane of the 
jarosite structure. 

In both the pure kagome and jarosite structures, the 
Dy are not completely determined by symmetry. How- 
ever, if we fix arbitrarily one of the Dy, then all the oth- 
ers are fixed by the three-fold rotation axis perpendicular 
to the kagome plane that passes through the center of the 
triangles of the kagome lattice. 

Moriya's rules constraint the Dy vectors with the help 
of the symmetry of the lattice, but they are not a proof 
of the existence of DMI in the kagome lattice (or the 
related jarosites), they just express the fact that if DMI 
are present, then the Dy will necessarily be restricted to 
some set of possible vectors. 



B. General considerations 

Taking into account the superexchange mecanism, the 

t 2 - 

isotropic exchange Jy is proportional to jj- (iy being 
the intersite hopping and U the on-site Coulomb repul- 
sion), while it was shown by Moriya that |Dy| is pro- 

Ai 2 - 

portional to (A being the spin-orbit coupling and 
A the crystal field splitting) and A y is proportional 



to 



it' 

"A T U' 



This last term, being one order of magnitude 
smaller than the DMI is often neglected. However, it 



was suggested by Shekhtman et al. Q that it plays an 
important role, as does the DMI, because they are re- 
sponsible for anisotropies respectively proportionnal to 
— and A which are of the same order of magnitude. This 
argument is based on the assumption that the isotropic 
exchange J, which is the dominant interaction favours 
collinear configurations, and is untrue in the case of the 
kagome lattice due to its frustration. In a collinear struc- 
ture DMI are in competition with isotropic exchange J 
resulting in a ^j- anisotropy, whereas anisotropic sym- 
metrical exchange defines easy axes or planes but is not 
contradictory with a collinear structure, resulting in a 
~A anisotropy. Considering different possible A ij, we 
find that these arguments generally do not hold in the 
case of a non-collinear structure such as those found on 
the kagome lattice. Rather, the anisotropies are respec- 
tively of the order of D and A, which is the reason why 
only DMI will be taken into account in this work. The 
anisotropic exchange A is considered to enter the prob- 
lem at the next order of perturbation theory. 

Shekhtman et al. also showed that under specific sym- 
metry conditions it was possible to map the total spin 
hamiltonian of Eq. (Q) onto an isotropic Heisenberg 
hamiltonian (only first term of Eq. (Q)) thus recover- 
ing a hidden rotational symmetry. These conditions are 
not fulfilled in the kagome lattice and this mapping is not 
possible. 



C. Microscopic derivation of DMI using Moriya's 
technique 

As mentionned above, while application of Moriya's 
rules restricts the D to some set of possible vectors, they 
are not proof of the existence of DMI. In this section we 
derive DMI using the method proposed by Moriya assum- 
ing some microscopic situation (orbitals and crystalline 
electric field) which is arbitrary, but related to the case of 
the jarosites. The method proposed by Moriya applies for 
localized magnetic electrons (insulators). In order for the 
derivation to be manageable, we assume that the ground 
state is not degenerate due to cristal electrical field (ex- 
cept Kramers spin degeneracy) and each magnetic site is 
occupied by one electron. In our case, we consider the 
t2 fl orbitals, whose degeneracy has been further lifted by 
some distortion of the octahedral environment. The e g 
orbitals are assumed to be at much higher energy level 
(or completely filled and at much lower energy level). 
The crystalline field scheme for the t2 fl is represented on 
Fig. [I]. The three states of each triangle are obtained 
from one another by a rotation of 4^ around the axis 
perpendicular to the kagome plane. 

The octahedral symmetry and the 3d orbitals are rele- 
vant for the jarosites (see Q), and the lifting of degeneracy 
of the t2 g orbitals is supported by the fact that the oxy- 
gen octahedra surrounding magnetic sites are distorted 
in the jarosites (see figures 2 and 3 of ref. |l5|). The local 
axes of the cristalline field on each site have also been 
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FIG. 1: Local axes of the crystalline electrical field considered 
in the microscopic derivation of the DMI. The axes are tilted 
by an angle a around the x axis toward the center of the 
triangle. 




tilted towards the center of the triangles of the kagome 
lattice by an angle a, in order to fit with the jarosites 
symmetry. Spin-orbit coupling is taken into account at 
first order in perturbation theory as it is expected to be 
much smaller than the crystalline electric field. Next, the 
on-site Coulomb repulsion is assumed to be much higher 
than the hopping term, and the intersite hopping is in- 
troduced to second order in perturbation. Doing so, we 
arrive at the following expression for the D12 vector 



D x 
D z 



V3D Y 
\V3 



256[/ 
' Ma) 
Ax 

AV3 



sin(2a)/i(a) 



cos 2 a 



32E/Ai 



— ((ddS) - (dda)) 
^2 

h ( a ) h( a ) 



where 



A(a) 
/2(a) 



4(ddTr) - 9(dda) sin 2 a 
-3(dd8)(l + 3 cos 2 a) 
4(cMtt)(3cos 2 a - 2) - 9(dda) sin 2 a 
+ (ddS)(l + 3 cos 2 a) 



(2) 

(3) 
(4) 

(5) 
(0) 



(dda), (ddir) and (ddS) are the transfer integrals as 
defined by Slater and Koster Q. 

The obtained D,y vectors are located in the plane per- 
pendicular to the (ij) bond and are perpendicular to the 
kagome plane if the octahedra are not tilted (a = 0) as 
they should according to Moriya's rules. 



III. CASE OF D _L KAGOME PLANE 

In this section the consequences of the DMI on the low 
temperature magnetic structure are explored in the case 
where D is perpendicular to the kagome plane. Both 
cases D z >0 and D z <0 (z is the axis perpendicular to 
the kagome plane and the convention is taken that the 
spins in the cross products appear rotating anticlockwise 
around the hexagons) are considered through mean field 
approximation and Monte Carlo simulations with classi- 
cal Heisenberg spins. 



FIG. 2: D _L kagome plane. The spins lie in the kagome 
plane and the sign of Dz selects the chirality. There is a 
global rotational degree of freedom around the z axis. 




FIG. 3: 1)1 kagome plane. The critical temperature is ~ D. 



Looking for q=0 structures, a restriction which will 
be justified by Monte Carlo simulations, it is found that 
within mean field approximation, the system undergoes 
a phase transition to a long range ordered state. In the 
ordered state, all the spins lie in the kagome plane, so one 
effect of the DMI is to act as an easy-plane anisotropy 
(as long as the first excitations are not taken into account 
e.g. looking only at the structure of the ground state). 
As it can be easily seen by expanding the cross product 
with only a z component for the D« , the hamiltonian is 
invariant under a global rotation of all the spins around 
the z axis. This degree of freedom is of course still present 
for the ground state. Since it is a q=0 structure, all 
triangles have the same magnetic stucture. Depending on 
the sign of D z , two chiralities are found and represented 
on Fig. § 

In order to study the behaviour of this system at finite 
temperatures, Monte Carlo simulations have been per- 
formed on finite size clusters with classical Heisenberg 
spins. In particular, the behavior of the critical temper- 
ature as a function of -j is plotted in Fig. ||. The critical 
temperature was taken as the maximum of the specific 
heat. 

The line D z =0 on Fig. || represents the well studied 
(classical) kagome lattice with antiferromagnetic near- 
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FIG. 4: In-plane D represented in the middle point between 
the magnetic sites and the associated magnetic structure. All 
the spins have a weak z component resulting in weak ferro- 
magnetism. There is no global rotational degeneracy. 

est neighbor exchange interactions. This system shows 
partial order: the coplanar states are assymptotically 
selected as temperature is lowered || (order by disor- 
der). However, as soon as D has a finite value, even 
much smaller than J, a phase transition occurs at fi- 
nite temperature. For this reason the low temperature 
magnetic structure will be governed by DMI. It must 
be emphasized that this argument would not be correct 
if DMI did not select a peculiar set of magnetic states, 
e.g. if the induced energy scale was renormalized to zero 
as it is the case for the exchange for the sole isotropic 
Heisenberg model. This is supported by the fact that 
the critical temperature is almost linear in D, and only 
weakly dependent on the strength of the antiferromag- 
netic exchange interactions (J). The physical reason for 
this is that the long range ordered states which are im- 
posed by DMI are part of the degenerate ground states 
of the system with only isotropic exchange (J) interac- 
tions, even after a partial lifting of degeneracy due to 
thermal fluctuations. Thus, DMI has here a first order 
effect on a degenerate ground state, while usually DMI 
acts as a small perturbation on antiferromagnetically or- 
dered ground state, leading to second order corrections. 



IV. CASE OF TILTED OXYGEN TETRAHEDRA 

In this section we study the general case where 
is in the plane perpendicular to the (ij) bond. In order 
to clarify the discussion, we start by looking at the sole 
effect of the in-plane component of D before turning to 
the general case. 



A. Case of D in the kagome plane : canted 
structure 

The D vectors as well as the ground state magnetic 
configuration obtained by mean field approximation and 
Monte Carlo simulations are represented on Fig. ||. 

The structure has some similarities with the D z >0 case 



for D_L kagome plane, and indeed it has the same chiral- 
ity, but there are also big differences. The spins do not 
lie in the kagome plane anymore, but they all have the 
same out-of-plane z component, giving rise to weak ferro- 
magnetism, perpendicular to the kagome plane. Chang- 
ing the sign of the in-plane component of D changes the 
sign of the z component of the spins as can be easily seen 
by expanding the cross product in the hamiltonian. A 
second difference is that there is no longer a global ro- 
tational degree of freedom and each spin points towards 
a fixed direction. In this case, DMI seems more to act 
like an easy-axis anisotropy. The ground state for the 
DMI alone does not belong to the ground state manifold 
for the exchange interactions, since the sum of the spins 
on each triangle is not zero and there is now a compe- 
tition between exchange and DMI. The consequence of 
this competition is that the angle between the spins and 
the kagome plane depends now on -j : this angle, as well 
as the weak ferromagnetic component, increases with -j. 

B. general case 

We now turn to the general case where Dy lies in the 
plane perpendicular to the (ij) bond. If D z >0 then the 
chirality which is selected by D 2 is the same as the one 
which comes from the in-plane component D p . Starting 
from the configuration for D z =0, the effect of D z on the 
lowest energy configuration will be very much the same 
as J: it will decrease the value of the angle between the 
spins and the kagome plane. The case D z <0 is more 
complicated : the chirality favoured by D 2 is not the 
same as the one which is selected if only the in-plane 
component of D is considered. The result is a competition 
between the in-plane and the out-of-plane components of 
D. When D is almost in the kagome plane with a small 
negative component along the z axis, then the canted 
structure is selected with a ferromagnetic moment which 
increases with increase of the z component. There is a 
critical value of jy 1 -, and if D 2 is negative enough, then 
the ground state is no longer the canted structure, but 
the planar structure on the right of Fig. as in the case 
of D _L kagome plane. This critical value depends on the 
strength of the exchange interactions J, since the latter 
favour the planar structure. 

It appears that the magnetic structures for D and - 
D are not always similar as in the case of D _L to the 
kagome plane. They can not be deduced easily from each 
other. This can also be seen on the hamiltonian : when 
replacing D by -D, there is no simple transformation 
of the Si which would leave the hamiltonian invariant. 
This comes from the fact that the lattice is not bipartite, 
since in a bipartite lattice, exchanging the two sublattices 
would transform Sj x Sj into — Si x Sj and leave the 
hamiltonian invariant. However, D p — * -D p corresponds 
to a trivial symmetry : changing simultaneously S z — > - 
S z will leave the hamiltonian invariant, and so the non- 
trivial part of the transformation D— >-D comes from the 



5 




0.7 
0.6 
0.5 
Zk 0.4 

J 

0.3 
0.2 
0.1 





-1.5 



S =0.5 



-0.5 0.5 1 



1.5 



6 



FIG. 6: Critical temperature as a function of the angle 
between D and the kagome plane, for a fixed value of 

T>=y/Dl+Di 



FIG. 5: Ground states obtained for different values of J 
(isotropic exchange), D p and D z (in plane and out of plane 
components of D). The intensity of the grey color represents 
77, the angle between the spins and the kagome plane. The 
grey on the left part corresponds to 77 = (spins in the kagome 
plane) while on the right part, the white corresponds to ??>0 
and the black to ?7<0. 



z component D z . The discrepancy between the structures 
for D and -D is illustrated on Fig. [| It represents the 
different ground states obtained for different values of the 
parameters D p , D z and J. 

This figure was obtained by comparing the energies of 
the two possible chiralities (a) and (b) , (see figure ||) on 
one triangle, and minimizing with respect to 77, the angle 
between the spins and the kagome plane and the angle tp 
between the projection of the spins on the kagome plane 
and their postion as represented on figure ^. The energies 
of the two triangles (a) and (b) are: 



When chirality (a) is selected, and D p 7^ 0, the spins 
have an out of plane component with the angle between 
the spins and the kagome plane being: 



2LL 

tan(2??) - 



VSJ + D Z 



(9) 



Ea 

N 



(1 - 3 cos (2??) ) - V3D Z COS 2 T] 



v3Dp sin(2?7) cos 
(7) 



— = -(1-3 cos (277)) 
N 2 v " 



V3D 



z COS 77 



(8) 



Several aspects already mentionned appear on these 
expressions. First, the sign of D z selects the chirality 
since it appears with a different sign in E Q and E;,. D p 
appears to favour chirality (a) and weak ferromagnetism 
perpendicular to the kagome plane (77 ^ 0). For chirality 
(a), if D p —Q, the spins lie in the kagome plane (77 = 0), 
and there is a global rotational degree of freedom since 
the energy would not depend on tp. If chirality (b) is 
selected, then 77 = always minimizes the energy (even 
for D p ^ 0) and so only in-plane structures appear with 
this chirality and the global rotationnal degree of freedom 
is always present. Furthermore, in this case the energy 
does not depend on D p . 



In this case, tp = (no global rotationnal degree of 
freedom) . 

Figure]^ shows the evolution of the critical temperature 
as a function of the angle 9 between D and the kagome 
plane. 

The critical temperature is lower when D is in the 
kagome plane (#=0) than when D is perpendicular (9 = 
±§). This is easily understood since for D in the kagome 
plane, DMI favours a canted structure, which is not a low- 
est energy configuration for the isotropic exchange part of 
the hamiltonian (J). For 9 = ±-|, DMI selects a coplanar 
structure, which also minimizes the isotropic exchange. 

The fact that the critical temperature is not symmetri- 
V cal with respect to 0=0 (Figure ||) also illustrates the dif- 
ferent behaviour for D and -D. Indeed, changing 9^—9 
corresponds to D z — > -D z which is the non-trivial part of 
the transformation -D, because D p — > -D p corre- 
sponds to the trivial symmetry S 2 — > -S 2 . The fact that 
the critical temperature is always different for symmet- 
rical values of 9 with respect to is due to the fact that 
D z <0 is in competition with the in-plane component of 
D since they tend to select different chiralities; on the 
contrary D z >0 and D p drive the system to the same chi- 
rality, hence a higher critical temperature is obtained for 
9>0 than for 0<O. 



V. MAGNETIC STRUCTURE OF FE AND CR 
JAROSITES 

We are interested here in the low temperature mag- 
netic structure of the Fe and Cr based jarosites, where 
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the magnetic ions (Fe 3+ and Cr 3+ ) form ka gom e planes 
which are stacked one on another || [ll], [l2|, [l| |l5) giv- 
ing rise to a three dimensional structure. However, the 
magnetic behavior of the jarosites is believed to be es- 
sentially that of a (two dimensionnal) kagome lattice, 
because the different kagome planes are fairly well sepa- 
rated by several non-magnetic atoms implying that the 
super-exchange interactions between neighboring planes 
are much weaker than between magnetic sites belonging 
to the same kagome plane. These magnetic sites are sur- 
rounded by octahedra of oxygen atoms responsible for 
the crystalline electrical field and involved in the super- 
exchange interactions. These octahedra are slightly dis- 
torted and tilted (see figures 2 and 3 of p5[ ) which is 
consistent with the microscopic electronic states we have 
considered for our microscopic derivation of DMI. 

A long range ordered state is observed experimentally 
in these compounds (^,11, [l2[ [l3[ |l5), and the low tem- 
perature magnetic structure has an ordering wave vector 
k=(0,0,0) or k=(0,0,|), depending on the diamagnetic 
atoms present in the compound. The third component 
of the k vector corresponds to the propagation of the 
magnetic structure from one kagome plane to another 
and is due to weak inter-plane super-exchange interac- 
tions which were not considered in this article since they 
play no role in the ordering of one kagome plane. 

Henceforth the ordering of only one kagome plane is 
considered. The different observed magnetic structures 
have a q=0 wave vector and always correspond to chi- 
rality (a). Both planar and weak- ferromagnetic (all the 
spins having the same z-component along the axis per- 
pendicular to the kagome plane) structures are observed 
experimentally and both are obtained here theoretically. 
The planar structure with chirality (b) obtained theoret- 
ically is not observed experimentally in these jarosites. 

The experimental structures can also be obtained in- 
troducing a single-ion anisotropy |ll| , [l3| , however this 



is not a relevant interaction at least in the case of Fe- 
jarosites. Indeed, the magnetic ions are Fe 3+ , and so 
the 3d orbitals are half filled, with five electrons coupled 
giving rise to a spin S=|. Thus, the charge distribution 
is spherical for this orbital occupancy, and consequently 
prevents the appearance of a spin anisotropy. Also for 
Cr 3+ ions, the 3 t2 ff orbitals are filled (S=§) and the 
single ion anisotropy should be small. 



VI. CONCLUSIONS 

It has been shown that DMI are allowed by the sym- 
metry of the kagome lattice and are relevant interactions 
whose effect is enhanced by the frustration of the pre- 
dominant antiferromagnetic isotropic exchange. A mi- 
croscopic derivation of these interactions for a schematic 
electronic configuration somewhat related to the jarosites 
structure was done. The magnetic properties due to 
DMI were studied through mean field approximation and 
Monte Carlo simulations. Contrarily to the predominant 
antiferromagnetic isotropic exchange, DMI can induce 
several long range ordered magnetic structures depend- 
ing on the different microscopic parameters J, D p and 
D z . The low temperature magnetic structures of the Fe 
and Cr based jarosites can be explained by the presence 
of DMI. This would be the first example of compounds 
where magnetic ordering is not due to exchange interac- 
tions but to DMI. 
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